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Abstract 



An extension of the Legendre transform to non-convex functions with 
vanishing Hessian as a mix of envelope and general solutions of the Clairaut 
equation is proposed. Applying this to systems with constraints, the pro- 

C|^ ' cedure of finding a Hamiltonian for a degenerate Lagrangian is just that of 

solving a corresponding Clairaut equation with a subsequent application 
of the proposed Legendre-Clairaut transformation. In this way the un- 
constrained version of Hamiltonian equations is obtained. The Legendre- 
Clairaut transformation presented is involutive. We demonstrate the ori- 
gin of the Dirac primary constraints, along with their explicit form, and 

N I this is done without using the Lagrange multiplier method. 

oo 

T^ ; 1 Introduction 



i>^ " Modern field theories are in fact degenerate dynamical systems whose key fea- 

''T^ . ture is the presence of constraints [Ij [2] . The most common way to deal with 

C^ '_ such systems is to use the Dirac approach "3; based on extending a phase space 

and constructing the so-called total Hamiltonian. In spite of its general success, 
e.g. in describing systems with gauge symmetries and gravity [HIS], the Dirac 
approach has limited applicability and some inner problems "B]. Therefore, it is 
worthwhile to reconsider several basic ideas of the Hamiltonian formalism per se 
starting from the Legendre transformation treated as a solution of the Clairaut 
equation also in the singular case. The Dirac approach is based on the follow- 
ing idea: in finding a Hamiltonian, to use the standard definition of momenta. 
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then perform the Legendre transformation and add Lagrange multipliers, sub- 
sequently to be removed by imposing some relations between constraints [3] . 

We revisit the procedure of finding a Hamiltonian for both regular and sin- 
gular cases. This procedure is reduced to that of solving the Clairaut partial 
differential equation. In the case of regular systems a Hamiltonian corresponds 
to its envelope solution [7] , while adding the general solution of Clairaut equa- 
tion leads to the total Hamiltonian of singular systems, and arbitrary constants 
correspond to the Lagrange multipliers within the Dirac approach [3]. Such 
solutions exist for smooth Lagrangian functions, while the standard Legendre 
transformation is applicable in the regular case only. To solve the Clairaut equa- 
tion in the singular case, we introduce a mixed envelope solution, which is an 
envelope solution in "regular" variables and a general solution in "nonregular" 
variables. We use the coordinate language which is convenient for making the 
basic idea transparent as well as for further applications. 

2 The Legendre transform and the Clairaut 
equation 

We recall the standard Legendre transform definition [7 . The Legendre trans- 
form of a convex functioru F : R" ^ M is a map £eg : F i — > G, where 
G (p) : R"* ^ R is another convex function (in the dual space) such that 

def 

G (p) = maxxG(p,x), and G(p,x) = p • x — -F(x). The maximum is at- 
tained, when 

9G(p,x) dF{x) 



dx d> 



0, (1) 



which, for a given p, can determine x unambiguously, but not vice versa. The 



de f 

convexity implies that the Hessian Hp (x) = det 



d^Fjx) 



is positive, which 



ensures that (globally) there is only one critical point, the maximum [7]. 

Now we reformulate the Legendre transform in terms of the differential equa- 
tion for G (p) as above. Let us suppose that ([T]) has a solution x = X (p) with 

X : R" ^ R", then 

G (p) "1=^ G (p, X (p)) = p . X (p) - F (X (p)) . (2) 

Since g 2 — ( a ) 0^2^ , the functions F and G belong to the same 
convexity class and in the case of convex functions one has rank q 2 



^We use vector notation 7 in which xSlRxRx...R = R", and x-y is the scalar product. 
Functions and variables are denoted by capitals and lowercase letters, respectively. A scalar 
differentiable function F : R" — > R is denoted by F (x). A vector dF/dx denotes a gradient 
of F whose entries are just partial derivatives dF/dx^, i = 1, . . . [i] = n, and for d^F/dx'dx^ 
we sometimes use the notation d'^F /dy? . Also, 9V/9x denotes the divergence of a vector 
function V : R" — + R" in a similar manner. 
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tlM 



rank — g^ 

us to exclude X (p) from ^ and obtain 



dGjp) Q 



Next we differentiate ^ and get q' = X (p) , which aUows 



which is just the Clairaut equation [5]. Obviously, ^ has solutions even in 
cases, where ([TJ cannot be resolved in x, and therefore the standard Legendre 
transform ^ does not exist (thus we add the superscript CI). We call a map 
£eg : F i — > G^^ defined by ([3]) a (generalized) Legendre-Clairaut transform 
with F (x) being any (smooth) function (not necessarily convex) . For a con- 
vex function -F'(x), if G (p) = Legp (p) is its standard Legendre transform, 
then G(p) satisfies ([3]), hence G(p) coincides with G*-^' (p) — Leg^'(p), the 
Legendre-Clairaut transform of F (x). Now we demonstrate the converse state- 
ment. Let us write a general solution of ([3]) as 

G^i„(p,c)=p.c-^(c), (4) 

where c G R^. The envelope solution can be obtained from the extremum 

condition — "T = p J^ = 0, which coincides with the condition ([1]) 

(this G^'„ actually coincides with G above). If Hp (c) > (i^ is convex), the 
extremum condition can be solved by c = X (p). Then the envelope solution of 
dSl) is 

Gfi (p) 'i^^ G% (p, X (p)) = p . X (p) - i^ (X (p)) i G (p) . (5) 

This means that, given a convex function F (x), if G^^ (p) = Leg'p'' (p), then 
G'-^'- (p) = Legp (p) = G(p). Also, only in this (convex) case both mappings 
£eg and £eg are involutive £,eg o £eg — id, £eg o £eg = id. Thus, the 
standard Legendre transforms in the class of convex functions are in 1-1 corre- 
spondence with the envelope solutions of the Clairaut equation ^ . This provides 
an exposition of the ordinary theory [71 [5] in a special way, which is convenient 
for our subsequent purposes that involve more general classes of functions re- 
lated to the Hamiltonian structure of constrained systems [21 IS] • 

Indeed, let us consider the case of a non-convex function F(x), when its 
Hessian Hp (x) vanishes. The standard Legendre trick does not work, because 
([1]) cannot be solved for x in this case [7]. On the other hand, the Clairaut 
equation Q assumes nothing but smoothness about -F(x). Therefore, we can 
forget the condition of its obtaining and start from the Clairaut equation itself, 
then try to find the solutions. In this way, we can extend the Legendre-Clairaut 



a^F(x) 



^ k < 



transform to the degenerate case Hp (x) = 0. Let rank 

and k is constant on the domain of -F(x). Without loss of generality, we can 
assume the indices are rearranged in such a way that a non-singular minor of 
rank k is in the upper left-hand corner. Then, we express the index i as a 
pair i — (11,^2), *i = 1, . . . /c, 12 = fc + 1, . . . n, and correspondingly any vector 
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variable is presented as x = (xi, X2), xi £ R'^, X2 G M""'^. We call the first and 
the second entry of x as regular and nonregular, respectively. In this notation, 
the Clairaut equation ^ becomes 

.ci, „^ „ 5GC'(pi,P2) _ 5GC^'(pi,P2) 



G^'(Pl,P2) = Pi ;"-"-^ +p2. 



9pi dp 



2 
^ 9G^'(pi,P2) dG^^Hpi^\ ff,. 

By analogy with Q , the general solution is 

<5^e«(Pl>P2,Cl,C2) =Pl •C1+P2-C2 -F(ci,C2), (7) 

where Ci G M*^, C2 G M"^'"'. Our intention now is to search for the envelope 
solution (O in regular variables only. But it is still the general solution ([4]) with 
arbitrary C2 with respect to non-regular variables. We call such a solution a 
mixed envelope solution of the Clairaut equation. Differentiate ([7]) in Ci to get 

c'G'^e«(Pl,P2,Cl,C2) aF(ci,C2) ^ ,„, 

d^, =P^ d^^^^ ^^) 

Since the sub-Hessian Hp ' (x) = det g^i'^^ does not vanish, we can 

resolve ([5]) with respect to Ci and obtain Ci = X(pi,C2). The subsequent 
substitution to ([7]) yields 

G^'(pi,P2,C2) - G^i„(pi,P2,X(pi,C2),C2) 

= Pi •X(pi,C2) -fp2 -02 -F(X(pi,C2),C2), (9) 

which can be treated as an explicit form of the Legendre-Clairaut transform 
for a non-convex function F (x) having degenerate Hessian matrix. Note that 
the relation between the Legendre transform and the parallel curves (which 
is connected with the above general solution) was considered a long time ago 

[smu]. 

Now we apply the Legendre-Clairaut transform to the Hamiltonian proce- 
dure for constrained systems with finite number of degrees of freedom (in the 
language of classical mechanics). This is sufficient for exploring the main idea 
which can be easily generalized to a field theory (e.g. using DeWitt's condensed 
notation |TT] . 

3 The Hamiltonian procedure and the Clairaut 
equation 

Let Q be a n-dimensional configuration space being a smooth manifold with 
local coordinates q = (q^,...^"') (all statements can be translated into the 
coordinate free language [IHIISj)- A Lagrangian on Q is a continuous function 
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C : TQ -^ M. that is smooth on the tangle bundle TQ\{0} which in local 
coordinates is determined bjo (q, v), where v (t) = ^^ ' are velocities. In this 
notation the Euler-Lagrangc equations of motion are 

where W(q),jj (v) = 9^>C(q) (v) /dv'^dv^ is the Hessian matrix and i^(q)i (v) = 

The standard Legendre transformation ^7 is a local mapping TQ -^ T*Q (the 
latter is the phase space which is (q, p) in local coordinates) or £eg : £ ^ H, 
where ?^ : r*Q — > R is a Hamiltonian. 

def II I 

First consider the regular case, when the Hessian Hc^^^ (v) = det || M/'(q)y (v) | 
is nonvanishing. Indeed, the observations of Section [2] imply that it is the Leg- 
endre transform in velocities (considering q as parameterf|j) such that 7i(q) (p) = 

def 

maxv7i(q) (p, v), where H(q) (p, v) = p • v — £(q) (v). The extremum occurs, 
when 

av P 9v ' ^ ' 

which can be resolved with respect to velocities v = V(q) (p) , since the Hessian 
is nonvanishing. Then in the regular case the Hamiltonian is 

^(q) (P) '= ^(q) (p. V(q) (p)) = P ■ V(q) (p) - ^(q) (V(q) (p)) . (12) 

Now we differentiate 7i(q) (p) and obtain 

9H(q) (p) IP 



dp 



V(q) (P) . (13) 



Because p^ holds for all q, p identically for a solution of (fTTj) , we are 
able to substitute V(q') (p) into (fT^ and obtain the Clairaut equation for the 
Hamiltonian as follows (cf. ([3])) 



(q)lPj=p-^g^— -Aq)(^9^ 



n.) (p) = p ■ ^^-^ - ^(q) ^^-^ • (14) 



We call this map a (generalized) Legendre- Clairaut transformation £eg : £ — *■ 
Ti,^'', because ([14]) has a solution also in the case of singular Lagrangians. In the 
regular case we follow the steps of the previous section in considering a general 
solution of pi]) 

Hfq)ge„(P,v)=p-v-/:(q)(v), (15) 



^As we consider time independent Lagrangians for conciseness, the time-dependent case 
can be treated similarly. 

^We write q-dependence as a subscript to single out coordinates as passive variables or 
parameters under the Legendre transformation for which v, p are the active variables. 
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where initially v = C(q) are constants with respect to the active variables v, 
p, i.e. arbitrary functions of q as parameters. The envelope solution of p3|) is 
subject to the extremum condition 

av av 

which coincides with (jll[) and determines additional dependence on the mo- 
menta when we resolve ([T6| (which is possible because He 7^ 0), and we denote 
this solution by v = ^(q) (p)- After substituting into (fT5|) we obtain the enve- 
lope solution of the Clairaut equation (fT4|) as 

W(q)e™ (P) = ^(q)sen (P> v) |v=V(q)(p) = P-V(q) {p)-£{q) (V(q) (p)) = H(q) (p) , 

(17) 
which coincides with the standard Legendre transformation (fT^ . as it should 
be in the regular case [8]. 

4 The Clairaut equation for constrained systems 

Consider a singular dynamical system for which the Hessian He vanishes. Direct 
application of the standard Legendre transformation is not possible now, because 
(fTTj) cannot be solved with respect to velocities v. But in the Clairaut equation 
p^ there are no restrictions on £(q) (v) except smoothness, and therefore, 
as in Section [2l we are able to consider the (generalized) Legendre-Clairaut 
transformation (fT4|) in the singular case He, (v) = as well. 

Let the rank of the Hessian matrix be less than the configuration space 
dimension ranfc WjqUj (v) = k < n and k is constant. We rearrange the 
indices i, j in such a way that a non-singular minor of rank k will be in the 
upper left-hand corner. Then, we express the index i as a pair i = (11,12), 
ii — 1, . . . fc, 12 = fc + 1, . . . n, and decompose sets of coordinates and momenta 
as q = (qi,q2) = ({?*' },({?'" }). P = (Pi:P2) = (fej , (fes}), calling the 
first and the second set as regular and nonregular coordinates/momenta, re- 



spectively. In this notation the Hessian matrix is W 



*^ ^ W W- 



/ \A/(11) \A/(12) \ 

( W(2i) W(22) ) ' '^^^^■'^ ^^"^ ^^ nonsingular det W(") 7^ 0, rank W(") = k. 



Then the Clairaut equation (|14p acquires the form 
H(q) (p) = Pi rr- +P2 jr- Aq) ^ ' ^ ' (1«) 



which we treat as a definition of Ti/l'-) (p) in the case of singular Lagrangians. 

We cannot derive this, as in the regular case, because there is no relation (|16p 
for nonregular variables. A general solution of this partial differential equation 
is 

^S)sen(P>Vl,V2) = Pi • Vi + P2 • V2 - ^(q) (vi , V2) , (19) 
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where Vi = Ci(q),V2 = C2(q) are arbitrary functions of the passive variables 
q. As compared to pB]) . we can find the envelope solution for the regular part 
only, i.e. we obtain the mixed envelope solution. The extremum condition for 
regular variables now is 






Pi 



^^(q) (VI.V2) 

9vi 



= 0, 



(20) 



which can be solved (due to detW^^^^ ^ 0) as vi — V(q) (pi,V2) |v2=C2(q) = 
V(q) (pi,C2(q)). Then we substitute this solution into (fTO)) and obtain the 
"mixed Hamiltonian" (or unconstrained Hamiltonian) in the form 



de/ ^Cl 



'^(<i)mixed{V:^2) - W(q)gen (PjVl,V2) |vi=V(q)(pi,V2) 
= Pi • V(q) (Pl, V2) + P2 • V2 - /:(q) (V(q) (pi, V2) , V2) , 



(21) 



where V2 = C2(q) remain arbitrary functions of the passive variables q. In 



this picture the relation (|20p is not a definition of the momenta, but rather a 
condition for the existence of the envelope solution for the regular part. A similar 
condition for the nonregular part does not exist. Therefore at this initial stage 
the nonregular momenta P2 have no connection with the Lagrangian (analogous 
to (|20p ). thus now a "true" (in the standard definition) phase space is formed 
by(qi,Pi)e(T*Q), only. 

Note that in [M] the passage from the Lagrangian /I^q) (v) to the general so- 
lution H/q^N (p, vi , V2) is called "a slow and careful Legendre transformation" , 

while the further passage to the mixed Hamiltonian 'H'P^-,^^^^^ (p, V2) is called "a 
reduction of the global Hamiltonian Morse family" . Also, the given Legendre- 
Clairaut transformation becomes exactly a generalized Legendre transformation 
as of [T5] . 

Now we consider the full differential of both sides of (pij) and use the ex- 
tremum condition (|20p . which gives 



^^fq)™^ed(P.V2) 



dqi 



^^(qW.ed(P'V2) 



9pl 

9^fq)rmxediP^^2) 
9p2 



V2-C2(q) 



■V2— C2(q) 



V2=C2 



9£(q)(Vi,V2) 



9qi 



'^l=V(q)(pi,C2(q)) 
V2=C2(q) 



+R!if (p,C2(q)) 
V(q) (pi,C2(q)) , 

C2(q), 



(22) 
(23) 

(24) 
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where 

<S(P'C2(.))= E ^i.yM-^^ (25) 



*(,) (p) = P2 ~ *(,) (pO , *(,) (Pi) = ^^'"^^ ^^''^'^ 



dvo 



, ^ (26) 

vi=V(,,(pi,C,(,)) ^ > 

V2=C2(q) 



and C2(q) are still arbitrary. Note that *(q) (pi) and therefore $(q) (p) have 
no dependence on the unsolved velocities V2 , because, if some of them appeared 
there, we could derive them from (P5|) . which contradicts the fact that the 
rank of the Hessian is fc. Then, using the Lagrange equations PH)) and qi = 
V(q) (PijC2(q)), q2 = ^1{n)^ wc obtaiu the "mixed Hamiltonian equations of 
motion" (or unconstrained Hamiltonian equations) 

d'Hrq)mixed (p) 



(q)™..e.lPJ ^ -pi+RW(p,q2), (27) 

qi, (28) 



9qi ^^---(q) 

d'H(q)mixed (p) 



9pi 

d'H{q)mixed (p) 

^^- - -.(q),p,,^..(^) 

d'H(^q)rnixed (p) 
9p2 

It can be shown that the system (l27|) - ((30)) leads to the equations of motion which 
are equivalent to the Lagrangian ones (jlOp . Observe that, if we use ([20|) also for 
nonregular velocities and momenta, we can present r.h.s. of the equation (|29p 



-*(q)(pi) + R[5(p,q2), (29) 

q2, (30) 



as — P2 + *(q) (p) + R/ { (p, q2)- Then (P7)) - ((5D|) becomes the standard system 
of Hamiltonian equations, when the following system of equations (generalized 
constraints) is valid 

R^ (P' ^2) = 0, *(q) (P) + RJS (p, q2) = (31) 

The generalized constraints (pij) are sufficient to introduce the standard Poisson 
brackets and the "correct" time evolution [7]. 
In the Dirac formalism [3J one imposes 

*(q) (p) = 0, *(q) (p) - 0, (32) 

which are the standard primary scleronomous constraints (which should be func- 
tionally independent, otherwise see [H]). Obviously then, the conditions in ((5^ 
are more restrictive than those in (|3ip . But, they lead to the "correct" phase 
space in regular and nonregular variables, in which for both sets of momenta 
the derivatives of the Lagrangian with respect to corresponding velocities, are 
initially treated as definitions [T2]. Indeed, this allows one to consider (q, p) 



^'^{q_)mixed (Pj 


<^2(q)) - 


0,(0) 


(Pl 


)+C2, 


<',(P.) 


= Pl 


•V(q)(p 


17^21 


;q)) 


-^M 


+ C2(q) • 


*(q) 


(Pl), 
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as points of the entire "true" phase space T*Q, while the Legendre transforma- 
tion becomes then a degenerate mapping with a kernel [12] . Then the "mixed 
Hamiltonian" ([2T|) can be presented as the sum of Ti, { (pi) and the linear 
combination of the primary constraints 

C2(q) • *(q) (P) , (33) 

V(q) (pi,C2(q)) ,C2(q)) 

(34) 

where hL! (pi) does not depend on p2 and C2(q) due to (pO)) and (|24)) . In this 
case ([33]) coincides with Dirac's total Hamiltonian [3] . 

If one does not impose the generalized constraints ([3T|) . then the equation 
([^^ has no time derivative of the nonregular momenta P2 on the right-hand side. 
Thus (|29|) is an algebraic equation for p2, and therefore there are no standard 
Hamiltonian equations for nonregular variables. Nevertheless, the system of the 
"mixed Hamiltonian equations of motion" describes the same constrained phys- 
ical system, as Lagrangian equations of motion (jlOp , independently of whether 
or not the constraints ([5^ are imposed. 

Finally, we can show that the Legendre-Clairaut transformation is involutive 
for constrained systems. Indeed, in full analogy with (|21[) . for a given "mixed 
Hamiltonian", satisfying ([27|) -(f30 l) . we can construct the corresponding Clairaut 
equation (similar to (fT8|) ). using the fact that the Hessian matrices for C and TC 
have the same rank (see Section [2]), from which we obtain its mixed envelope 
solution. This can then be called a "mixed Lagrangian" (or unconstrained 
Lagrangian) 

^{^)mixed (V, P2) = Vl ■P(q) (Vl, P2) +V2 ■P2 -Hfq')„i^,<i (P(q) (vi, P2) , P2, V2) , 

(35) 

where P(q) (vi,p2) = pi is a solution of the condition vi — — '''''"g''^'' — ^ — 
similar to ^ , and P2 are arbitrary functions (unsolved momenta) . If the func- 
tions P(q) (vi,P2) and V(q) (pi,V2) are mutually inverse in regular variables, 
i.e. P(q) (V(q) (pi,V2),P2) =Pi, V(q) (P(q) (vi , pa) , V2) = Vl, then it can be 
verified that /^mixed = ^^ which proves involutivity. 

5 Conclusion 

To summarize: the concise treatment of constrained systems using the Clairaut 
partial differential equation presented here gives a different explanation of the 
Dirac primary constraints from the Lagrange multiplier method, as well as a 
different understanding of the nature of primary constraints. In some cases this 
can lead to possible generalizations. Moreover, it can be applied in the cases 
where the standard constraint methods do not work or are too cumbersome. 
Examples, details and corresponding Hamiltonian formalism will appear in the 
forthcoming paper. 
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